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Abstract Three-dimensional central symmetric bodies different from spheres 
that can float in all orientations are considered. For relative density p — \ there 
are solutions, if holes in the body are allowed. For p ^ \ the body is deformed 
from a sphere. A set of nonlinear shape-equations determines the shape in 
lowest order in the deformation. It is shown that a large number of solutions 
exists. An expansion scheme is given, which allows a formal expansion in the 
deformation to arbitrary order under the assumption that apart from x — 0, ±1 
there is no x, which obeys Pp^2{x) = for two different integer ps, where P are 
Legendre functions. 

1 Introduction and Summary 

A long standing problem asked by Stanislaw Ulam in the Scottish Book [T] 
(problem 19) is, whether a sphere is the only solid of uniform density which will 
float in water in any position. Such a solid is called a floating body of equilibrium. 
It will be in indifferent equilibrium in all orientations. 

The simpler, two-dimensional, problem to find non circular cross-sections of 
a long cylindrical log which floats without tending to rotate (the axis of the log 
is assumed to be parallel to the water surface.) was solved for relative density 
Pd = 1/2 in 1938 by Auerbach [2 and for densities pd 7^ 1/2 by the present 
author [3 HIE]. 

Here we start to investigate the problem for three dimensional systems. In 
section 2 the basic equations for a partly immersed body in indifferent equi- 
librium are derived. In the following section the theorem that a star-shaped 
inversion-symmetric body in arbitrary dimension d and density 1/2 is a sphere, 
is reconsidered. One easily sees that a floating body of equilibrium, which is 
not star-shaped can have a shape different from a ball. A class of such bodies 
is explicitly given. 

The following sections, which are concerned with star-shaped floating bodies 
of equilibrium with arbitrary relative densities, can be read without having 



read section 3. In section 4 and 5 we start considering the expansion of the 
surface of a floating body of equihbrium around the sphere. To perform the 
expansion the distance of the surface is measured from the center of gravity. 
It is expanded in spherical harmonics, which are defined as the eigenf unctions 
of the Laplacean operator on the unit sphere with eigenvalues —1(1 + 1). For 
given I the eigenfunctions are linear combinations of 2/ + 1 linearly independent 
functions. 

In section 6 the equations are considered 
for the deformation in first order. If the angle 
6q is a zero of the associated Legendre function 
Pp,_2(cos 6'o) ~ 0, then the spherical harmon- 
ics with I = p contribute to the deformation. 
This equation is fulfilled at 6*0 = tt/2 for all 
odd p, which corresponds to pd = 1/2. Thus 
one has to expect that at this special density 
there is a large set of solutions. 
In the following we assume that with the ex- 
ception of cos 00 = 0,±1 there is no Oq, which Figure 1: Sphere of radius 
solves Pp__2(cos6'o) = for two different inte- ro and water-line for the defi- 
ger p. nition of h, p, and ^o- 

In section 7 the equations in second order in the deformation are derived. 
One finds that only spherical harmonics for the deformation are allowed whose 
square projected onto the space of harmonics with the same p is proportional 
to the original deformation. This is the contents of the shape equations 
This projection vanishes for odd p. Thus one has to distinguish between the 
problem for odd and even p. In the following only the case of even p, that is for 
central symmetric bodies, is pursued. A large number of solutions for the shape 
equations is found in section 8 assuming invariance under various subgroups 
of the orthogonal group 0(3). By considering only these special groups only 
solutions with mirror symmetry are found. It is not clear whether there are 
shapes with inversion symmetry, but without mirror symmetry. 

In section 9 it is shown how contributions in higher order can be obtained, 
which lead to a formal expansion for the shape. 




2 General Considerations 

Potential energy Denoting the volume of the body by V, the volume above 
the water by Vs,, that below the water by Vb, one obtains (Archimedes' law) 

- (1 - Pd)V, Vb = pdV, (1) 

where pd is the relative density of the body with respect to the liquid. Denote 
the total mass of the body by m, the masses above/below the water-line by 
ma,bi the center of mass above/below the water-line by Ca,b, and the distance 
of Ca,b from the water-line by (ia,b- Then the potential energy V of the system 
is given by 

V = rna.gda + (m - mb)gdb = m(l - Pd)g{ds. + db)- (2) 

Thus the difference in height between the two centers of mass, d = da + db, is 
constant, since it has to be independent of the orientation. This does not imply, 
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that da and db are separately constant. Moreover the line CaCb connecting the 
two centers of mass has to be perpendicular to the water-level. Placing the 
center of mass of the body in the origin and denoting the coordinates of Ca,b 
by (0, 0, Za,b) one obtains 

= Pdrf, Zb = -(1 - Pd)rf- (3) 

Thus the loci of the centers of gravity lie on spheres. 

The reverse is also true: If Za is independent of the orientation, then also 
Zb has this property. Then Za — Zb = da + db is constant and thus its potential 
energy. We will use this property to determine bodies which can float in all 
orientations. 

One can conclude: If a body has the property that a plane in arbitrary 
orientation cutting through it, so that the two volumina I4^b are constant and 
the centers of gravity of these volumina lie on spheres, then it is a floating body 
of equilibrium, since the body assumes in all orientations the same potential 
energy. 

These properties and another one we will derive below were derived by Pierre 
Bouguer and Charles Dupin a long time ago. Auerbach writes in his 1938 
paper [2]: II resulte aisement des theoremes classiques de Bouguer et Dupin que 
la condition necessaire et suffisante pour qu'un corps soit une solution de ce 
probleme est que la surface de centres de carene soit une sphere. De plus, on 
peut affirmer que dans ce cas V ellipse centrale d'inertie de la flottaisson est 
un cercle dont le rayon est le meme pour toute position d'equilibre. Indeed 
Pierre Bouguer (1698-1758) and Charles Dupin (1784-1873) wrote books[i [7] 
on the hydrostatics of ships. These conditions are also found in the textbooks 
by Appell[8] and Webster [9]. Compare also the article by Gilbert [lO]. 



Moment of inertia Rotate the body by an infinitesimal angle. Since the 
volume above and below the water has to be conserved the rotation is around 
the center of gravity M of the water-plane area F (intersection of the plane of 
the water-surface with the body). The center of gravity M of the water-plane 
area is given by 

XM = dFx, yu = — dFy. (4) 



F J F . 

An infinitesimal rotation S(j)x,y will bring a wedge of thickness 

({x,y) ^ d(j)x{x - xm) + S(f>y{y - yu) (5) 

above the waterline, if positive; if negative, its modulus describes the thickness 
of a wedge disappearing below the waterline. There are two contributions to 
the shift of the centers of gravity in x and y direction, 

5xa = -zj(j)x + ^ I dFxC(x,y), 5y^ = ~zj(py + ^ j dF yC{x,y),{Q) 



Sxh = -ZhStj): 



- -^y" dFxC{x,y), Syh = ~ZbS(f>y ~ ^ J dFy({x,y)l7) 

where the first term is due to the rotation of the centers of gravity, and the 
second one comes from the appearance and disappearance of the wedges. The 
requirement Sx^ — 6xb, Sy^ — 6yb yields for arbitrary S(l)x,y the conditions 

hj=6,jl, I^pd{l-pd)Vd (8) 
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for the moments of inertia of the water-plane area 



hx = jdF{x-XMf, (9) 
Ixy = J dF {x - XM)iy ~ Vm), (10) 
lyy = fdFiy-yuf. (H) 



The moment / is independent of the orientation. 

If ^ is not fulfilled, but both eigenvalues of the matrix 

/ = f j"" j"^ ) (12) 

are larger than /, then the body is in stable equilibrium. If both are less, 
then the body is in unstable equilibrium, and if one is larger, one less than /, 
then it is in saddle-point equilibrium. Here we are interested in an indifferent 
equilibrium, for which equality ([S]) holds. The corresponding condition yields 
in two dimensions that the length of the waterhne has to be constant. [3] 



3 Central Symmetric Case pd= 1/2 

In this section central symmetric bodies with relative density pd — 1/2 are 
considered. Central symmetry, which is also called inversion symmetry, means: 
If the point at r belongs to the body, then also the point — r belongs to it, 
where the center has been placed at the origin. The plane of the water-plane 
area goes through the origin in all orientations, since it cuts the body into two 
equal halves. 

The following theorem due to Schneider [TT], fT^ and to Falconer|13j. also 
referred to by Hensley in the Scottish book[l] holds: For arbitrary dimension d 
and density 1/2, if the body is star-shaped, symmetric, bounded and measurable, 
then it differs from a ball by a set of measure 0. It follows from theorem 1.4 
of Schneider [T^] (similarly corollary 3.1 in [TT]): 'Let fid be the unit sphere 
|u| = 1, and • the inner product. If $ is an even real-valued, countably additive 
set function on fid satisfying /^^ |u • v|d$(u) = for each v G fid, then <I> — 0.' 
The theorem will be reproved in the following assuming that the distance of 
the surface from the center of the body is a continuous function of u. This is 
done first for d = 3 dimensions and then generalized to dimensions d > 3. The 
premise that the body is star-shaped is important. It is shown in subsection 
p.3p that there are non-spherical floating bodies of equilibrium, if one allows 
holes to be drilled into the body. 

3.1 Star-shaped body in three dimensions 

If the body is star-shaped, i.e. there exists a point A such that for each point 
P in the body the segment AP lies in the body. Since the set of these points 
A (called kernel) forms a convex set and since it is central symmetric, too, the 
origin is such a point A. We denote the extension of the body from the origin 
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in direction of the unit vector u by r{u). Let flu be the unit sphere and v the 
normal on the water surface, then 

Ivd = J dVdf]„r|u • v| = J y dr!„r4(u)|u • v| (13) 

has to be independent of v. One expands |u • v| in Legendre polynomials 

I I n / N . N„+i(2n+l/2)r(n- 1/2) , ^, 

|u.vH;^cA(u.v), c„^(-)"+^ ^ 20F(n + l)! • ^''^ 



n=Q 



Note that all coefficients c„ differ from 0, c„ 7^ 0. Using the addition theorem 

47r 

n.m n.m \ ^ J : V ^/ 

Zn + I ■'^ — ' 

m 

where the l2ri.m are An + 1 orthonormalized (real) spherical harmonics, one 
obtains 

^ 2n + 1 J 

n,m 

Vd is only independent of v, if all coefficients c„ except cq vanish. Since the 
functions l2n,m form a complete orthogonal set of central symmetric functions, 
this implies that r^(u) has to be independent of u. 

3.2 Star-shaped body in higher dimensions 

The proof is similar to that in d dimensions, and r'*/4 have to be replaced by 
r'^ and r'^+^/((i + 1), resp. Instead of Legendre polynomials one has to expand 
in ultraspherical (Gegenbauer) polynomials G^f^^ ^''(u • v), 

oo 

|u.v| = ^c„Gr/^-)(u.v), 

11=0 

(-)"+H2n + d/2 - l)r(n - l/2)r(d/2 - 1) 
" " 27rr(n + d/2 + 1/2) ' ^ ' 

Repeated application of Gegenbauer's addition theorem (see appendix[X| yields 

2Tr^ 

(n + d/2- l)r(d/2- 1) 

_e functions Yn 
obeying 

A„y„W (u) - ^n{n + d- 2)Yj^%{u), (19) 

with the Laplacean operator A„ on the sphere fid- Since all c„ 7^ the same 
argument applies that r''+^(u) has to be independent of u. 



G(d/2-i) (u . v) = — V r^'^) (u)y('') (v) (18) 

lu vj („ + rf/2_i)r(d/2-l) ^ "'"^^ ^ ' 

where the functions Yn'^m form a complete set of orthonormal functions on fid 
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3.3 Bodies with holes 



If the requirement that the body is star-shaped is omitted, but holes in the body 

arc allowed, then the condition that the body K can float in all orientations 
requires that the integral /q drx{r,u)r'^ does not depend on the direction u, 
where the indicator function 

is introduced. Thus one may drill holes into the body and add the corresponding 
amount of material at the outside creating a bulge. 

Figure 2: Non-spherical central 
symmetric floating body of equi- 
_Q ^ librium with two holes described 

" ' in this section at = 1/2 pro- 

rh^ro=0-l jected onto the (x,y)-, (x,z)-, and 

(y,z)-plane. 

The inner part of the cylindrical 
holes end at the sphere, the outer 
, part at the bulge. 

y \ 




z 




As an example consider two circular cylindrical holes of radius rh around 
the axes x = ±b, y = drilled into a sphere of radius Ro. Using spherical 
coordinates 

X = rsinO cos cp, y = rsinO sin 0, z = r cos 6 (21) 
the bounds of the hole is given by 



b cos ± A / r? — 6^ sin 6 

r = R± = . (22) 

sm 6* 



One obtains the indicator function 



X{r,tf,9) I eiRo-r) otherwise 

(23) 

with the Hcaviside step function 9, which equals 1 for positive argument and 
for negative argument. 9{0) may be chosen 1 (0) for a closed (an open) body. 
R' is given by 

R' = ^R^ + RX-Rt. (24) 

Obviously there are bodies with several holes and the holes necid not be c;ircular 
cylinders. Thus there is a large variety of bodies with central symmetry at 
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density pd ~ 1/2, if one does not require star-shape. The body described here 
has the property that any fluid above the center can flow off and any air below 
the center can climb up up. In two dimensions a body with this property is 
not possible, since two holes through the body would cut it into three pieces. 
Otherwise, however, fluid above the center cannot always outflow nor air below 
the center cannot always bubble up. 



4 Parametrization of the Surface 

In the following bodies of general densities are considered, however, star-shape 
is assumed. Thus in spherical coordinates the radius r of the surface of the body 
is a unique function of the angles 0k, 0k in the body- fixed system. 

It is expanded in the set of spherical harmonics Ci^raiQ^)^™'^^ (Racah uses 
this normalization for C, however here C contains only the ^-dependence) , 

r = 5I«'.™^'^™(^K)e""*-, (25) 
which are related to the conventional spherical harmonics by 



Fz,„(0k,0k) = ^^^G,™(^^K)e™^«. (26) 

In contrast to section[3]and appendix [XI where the spherical harmonics were as- 
sumed to be real, here the spherical harmonics conventionally applied in physics 
are used dH [121 [H]. In particular one has 

Coo = l, Cio((?) =cos(0), Ci,±i(0) =T^sin0. (27) 

and in general 



G„.(0) = (-)"J^5— ^P/,„(cos0) (28) 
y (/ -I- to)! 

holds with the Legendre functions defined by 

1 Al+m 

P,„(.) = (l~.^)W2___(,2_iy. (29) 
The following relation holds 

Ci-^[d) = {-TCUB). (30) 

The body is rotated by means of Euler angles (a, /?) into a position where the 
water-surface is parallel to the x, j/-plane. A third rotation around the axis per- 
pendicular to the water-surface is not necessary. The functions C;„i' (f?K)e"" 
transform according to 

Q™K^K)e""'^« =^i?Ln'(«,/3,O)Q™(0)e™^ (31) 

771 

from which we obtain the radius 

r(0,</);/3,a) = ^ G™(0)i^^„, (a, /?, <^)a,,™. (32) 

l.Tfi.ra' 
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with the Wigner rotation matrix elements [THl HH [H] 

Dl^, (a, /3, 0) = e'™^d5,_, (/3)e™'". (33) 

The height of the water level above the center of mass of the body be h{(3, a). 
Then the intersection of the water level with the surface of the body is given by 

= cos(e)r(e,(/);/3,a). (34) 

This determines /3, a) of the waterline on the surface of the body for given 
orientation (a,/?). At water-level the waterline is given in polar coordinates 
(P, 0) by 

p(0,/3,a) =sin(e)r(e,0;/3,a). (35) 
The total volume V and the volume T4 above the water surface are given by 

V = \ j ^'^ der^sinO, (36) 

1 f /■®('^) h r 

K = 3 7 ^"^/^ dOr^sine-- J d(l)p^. (37) 

The volume 14,, which is a segment, is divided into a sector (first integral) 
and a cone (second integral). Obviously one has Vh = V — Va- Similarly one 
obtains the coordinates z of the centers of gravity by means of Zc = Vz^ and 
^a,b = K,bZa,b from 

= de'r'* sin 6* cos 6*, (38) 

= 4^ dOr'^ sine cose ~ — j d(j)p^. (39) 

The integrals p6ll39p can often be expressed in the form 

lm7n' Im' 

= 2^^G™-(/3)e™'"/z™.. (40) 

Im' 

Since the integrals must not depend on of a, /? all /;„i' have to vanish with the 
exception of /oo- 

4.1 The Sphere 

For a sphere of radius rg with the waterline aX Q — e^ one obtains (see figure 1) 

/i = rocos0o7 P^J'osin^o, (41) 

independent of the orientation. Since the cosine and sine of e^ is often needed, 
the abbreviations 

Co = cos0o, So = sin^^o (42) 
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are used. The volumes are 



3 

and the ^-coordinates 



K - -ri^(2-3co + c3) = -r3(l-co)2(2 + Co) 

= ^{2rl-irlh + h^), (43) 
Vy, = |r3(2 + 3co-c3) = |r3(l + co)2(2-Co) 

^{2rl + irlh-h^), (44) 



3ro(l + co)2 3ro (1 - Cq)^ ^^^^ 



One easily calculates 



4 2 + Co ' 4 2-Co 



3ro r 4 



d=Za-Zb = -— r, I=-rP (46) 

(2 + Co)(2-Co) 



and verifies eq. 



5 Expansion 

Starting out from the sphere the body will be deformed and an expansion in the 
deformation parameter e will be performed. That is, a and similarly the other 
quantities are expanded in powers of e, 

az.m = ^e"a„;i,m- (47) 

n 

Starting point of the expansion is the sphere of radius tq, 

ao;i,m = ro5ifl5mfi. (48) 

Similarly r is expanded, 

r(0,0;/3,a) =^e"r„(0,0;/3,a) (49) 

n 

with 

r„(0,0;/3,a) = ^ G™(0)Z?J„„, (a, A 0)a„;,„, (50) 

Since has to be evaluated around = Oq^ the expansion 

r„(0,0;/3,a) -;^((?-0o)V„fe(</),Aa) (51) 

A; 

with 



Inira' 



^mm'(a:P;0)a«;;m' (52) 
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is used, ro does not depend on 6, (j), (3, a. Similarly Q{^), h and p are expanded 
in powers of e, 

e(</.;/3,a) = 5^e"e„(,^;/3,a), (53) 

n 

HP, a) = ^e"/i„(/3,a), (54) 

n 

p{<j>,P,a) = 5^e>„(<^,/3,a). (55) 

n 

Qo — (^0, ho — fQCo, po — roso do not depend on the angles, since the body is a 
sphere in lowest (zeroth) order. 

5.1 Evaluation of the Integrals 

In order to evaluate the integrals for the part of the body above the water-level 
I divide the integral over 9 into one which runs up to and the one from 6 = Oq 
to 0((/'). Let us first consider the integral up to ^o- For the volume we expand 
in powers of e and in spherical harmonics 



^3, 



V3 = ^ e"/„ = e"a,n(^)e™*/„;;™. (56) 

n nlm 

Due to the (^^-integration only terms with m = will contribute. With ( = cos 9 
one obtains 

r d9Cio{e)sme= [ dC PiiO (57) 

Jo Jco 



yW^P;,_l(C) =-SoPi,-i(co). (58) 

COS 



Thus the first contribution 14*^^ reads 



1) = y" d</. ^ " sin 0/„ = 27r ^(1 - co)/„;oo - so ^ /n;(oPi,-i (co) j . (59) 

The term for Z = is written separately, since Po,-i is not defined. The total 
volume yields 

Va = 47r/„.oo. (60) 
Next the corresponding integral for Z is considered. The expansion 

rV4 = J2^"fn = E^"^'-W«""'^/n;<-- (61) 

n nlm 

yields 

[ ° d9Cio{9)sm9cos9 = [ dCC^i(C) (62) 

Jo Jco 

= cVi^Pi,-ii0' - [' ^/l^Pl,-l{0 

CO Jen 



= (CVW^^;,-i(C)-(i-C')i^/,-2(0) (63) 

Co 

= -SoCoP;,-l(co)+sgPi,-2(co). (64) 
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Thus the first contribution to ziV reads 



= I dckj'j d9f^ sin ^ cose = 27r Q(l - c^f^.^o + ^(1 - cg)/^;io 
+ so^/^.;o(soP(,-2(co)-coP,,-i(co))) . (65) 

l>2 J 

For the complete volume one obtains 

Zen = f /^;10- (66) 

We require that the origin coincides with the center of mass of the complete 
body. This yields /^-i o = 0- This can always be done by an appropriate choice 
of the an;i,m- 

Since Q{(f)) — 9o = 0{e) the integral over the interval = 6o...Q{(p) can be 
expanded in powers of e. Thus the second contributions to Vg, and Z^, yield 

V^n^ = J d<^<?„ (67) 

with 



dOr^sine- -p"^. (68) 



6' 

Similarly the second contribution to is obtained from 

Zi^ = I d.^5^ (69) 

with 



1 r^W 



d6lr* sin 6 cos 6 --h^p\ (70) 



It is useful to perform the expansion procedure explicitly in zeroth, first and 
second order. Then it becomes apparent how to continue to higher orders. 

5.2 Zeroth Order: The Sphere again 

In zeroth order one obtains 

/o = \rl = \rt, (71) 

90 = -^rlcosl, = -^r^cgsg. (72) 
One obtains from these expressions 

V^^^ = 2^(l-co)/o;oo = Yr^(l-co), (73) 
V^' = -^rlcsl (74) 



11 



which yields the volume in agreement with (j43p . and 

Z^^^ = 4l-c^)/o^^oo = J(l-cg>o, (75) 
^1? = -^il-cDclrt, (76) 
which yields Za in eq. (|45p . 



6 First order in the deformation 

The deformation of the sphere in first order in e is considered in this section. 
One obtains from eqs. (|34|35p 

©1 = , Pi = • (77) 

foSo So 

Following the procedure in subsection lS.ll one obtains fi = r^ri or more explic- 
itly 

h=J2 /i;i™Q„(0)e™*, = rlJ2 ^™U'(«- O)"!;'™' (78) 

and thus the volume 

Vi = Anruoorl (79) 
We choose to keep the volume unchanged in this order, which implies 

ri-ao = ai;00 = 0. (80) 
Next one obtains gi = — ^ft-irpSo and thus 

Vai = -27rsoro ^ ri;;oPi -i(co) - nhirlsl. (81) 
i>i 

Since T4i is constant, one obtains 

hi = hi-oo > ^ n-ioPi-i (co) (82) 

So 

with the a,/? independent contribution /ii;oo- 

Next one determines Z. With /f = rpri one obtains 

Zci = Y^o^iao- (83) 

Since the center of gravity stays at the origin, one requires ri-iQ = 0, which is 
equivalent to a^im = 0. With gf — —^hirgCoSQ one obtains 



Zai =7rSorj] 2^ri;;o(soP; -2(Co) - CqP; -i(Co)) - /iiCqSq 

V l>2 J 

Inserting eq. into the expression for Zai yields 



(84) 



Za.1 = TTSgro I 2^ri;/o-Pi,-2(co) - /ii;OoCo ) (85) 



l>2 
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Zai must be constant. We note that 

ri-iQ = ^ Dg^ia, (3, Q)ai.im- (86) 

m 

In order that Zai is constant aU coefficients ?'i:;oP;^_2(co) except for I — have 
to vanish. Therefore for any given / either all coefficients ai-im or P;,-2(co) have 
to vanish. We note that P;__2(±l) ~ for all I. However this corresponds to the 
limits pd = and pd = 1, in which case only a sphere is possible [17|. Further 
one has i-;,-2(0) = for all odd I. This corresponds to pd = 1/2. It may well be 
that similarly to the two-dimensional case one has a large variety of solutions 
for this special density. This case will not be considered further in this paper. 

In the following only cases are considered, for which Cq is different from 
and ±1. We will assume a cos 6*0 = Co for which -Pp. -2 (co) = and assume that 
for this Cq there is only one solution p. (Numerical calculation up to Z = 100 
shows that apart from Cq — 0, ±1 there is no Cq, which for two different ^s yields 
P/ _2(co) = 0.) Thus all ai-im — with the exception oil = p and one obtains 

ri =Y.'^pm{0)Dl^,{a,P,4>)a^,pm', (87) 

m 

where the amplitudes ai-pm' are as yet undetermined. Similarly /ii;oo is unde- 
termined. Note that the zeroes of Pp.2{x) and Pp,-2 are identical. 



7 Second Order in the deformation 

One obtains from eqs. ([Ml ISS]) in second order in e 

62 = — 4-3(-2so/i2ro + 2coSor2or-o - 2coSor'ii/ii 

-H2c^Soriirio - cohj - Co(2 - c^)r?o + 2hirio), (88) 

P2 = 7^ — 3-(-2coSo/i2?'o + 2sor2o?'o - 2soriift,i 
zroSQ 

H-2soCoriirio - hj - cgr^g + 2co/iirio). (89) 

These quantities depend on 0, a and p. Quantities on the right hand side 
depending on 6 have to be evaluated a,t 6 — do. The deviation of 9 from 6*0 is 
taken into account by the derivatives ru with respect to 0, compare eqs. ([51] 
[SI]). For products 6(1^6(2) of 

b^'^ =J2b'iilCi„^i9)e'"^'^ (90) 

Im 

we use the notation 

Mi)6(^)=5](6(i)6(^));i™G™(^)e'™^ (91) 

1 771 

with 

{b^'^b^%im^ J2 {I0;li0j20){lm;hm,j2m2)b[ll/^ll^, (92) 

h ,mi ,(2 ,m2 
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where the {Im;limi,l2m2) are Clebsch-Gordan coefficients. Use is made of 

=^(/™;?i™i,;2m2)(/O;ZiO,;2O)G,™(0). (93) 

Ira 

Then one obtains from 

/2 = rorl+rlr2. (94) 
/I = rlr^ + lrlrl (95) 

the expressions 

/2;im = rQ'^D''^^^,(a,l3,Q){{a\).im' +roa2-im'), (96) 

m' 

/I;/™ = '^oEA2™'(«,/5,0)(|(a?);i™,+roa2;,,„0- (97) 



For the complete volume one obtains 

2 , _[Q 

2p 

47r 

T 



V2 = 47r/2;oo = 47r(roa2;oo + '^{-)"^ ai;pmai-p-,n) , (98) 

/l;i,o- (99) 



Since the center of mass should stay at the origin and since (aD-im' vanishes for 
odd I one has (of );im = and requires a2:im = 0. 
Then the first contributions to Va2 and Za2 yield 

K2^ = 2^ |^(l-co)/2;oo-so5Z/2;ioPi,-i(co)j , (100) 

^(l-cg)/|;00 + ^i 



= 2^( ^(l-cg)/|^oo + ^(l-c;^)/2,io 



+ SoY^fhoMM^o) - CoPl,-iico)) I (101) 

;>2 

with eqs. (PI [^1571). 

The second contributions are obtained from 

92 = '^{-slh2ro + corl-2rihi+Cohl), (102) 
92 = ^rl{-2coslh2ro + 2clrl-icorihi + {3cl~l)hl). (103) 

together with eq. ([H71 15^ . The condition that V"a2 is independent of a, /3 yields 
an equation for Thus h2 is determined up to some additive a, /3-independent 
contribution /i2;oo- Instead of considering Za2 it is more practical to consider 
the difference Za2 — coroVa2 

'(2) „__T/(2) _ _ „ , 



- CoroK?' = const - 2nrlPl_M ^ i?L(a?);;^ 

+ 27rrQft.i;ooSo/'p-i(co)E£'om"i;p™- (104) 
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Then one obtains 



^a2 - CoroT42 - const (105) 

+ A2ih + ^3;(ai);;m) 

l>2,m 



with 



An = s2p,,_2(co), (106) 
A21 = so6pjPp^^i{co), (107) 

A31 = ^s2p,^„2(co)-^soCoP,,-i(co)-Pp%i(co). (108) 

From these expressions a2am is determined for Z 7^ 0, l,p 

3 1 

roslPl-2{co)a2:l7n = (- -SpP; -2 (cq) + ^CoSo^i,-! (Cq) + -Pp,- 1 (Cq ) ) («? ):im ■ 

For / = p one obtains the equation 



(109) 



So/ii;Ooai;pm = (^CqSo + Pp-i{cQ)){al).p„i. (110) 

For odd p the only non-trivial solution can be obtained for /ii;oo = Oi since the 
right-hand side vanishes for odd p. For even p a set of quadratic equations 

7ai;pm = ^ {pm;pmi,pm2)ai.pmiai.pm2 (HI) 

mi ,7712 

is obtained with the constant 

^ so/j-i;oo (112) 

(icoso + Pp „i(co))(pO;pO,pO) ' 

I call the equations (|llip shape-equations. They determine the possible shapes 
of the body in first order in e. Descriptively it means that the projection of the 
square of the deformation ai onto the harmonics with I = p is proportional to 
the deformation. I expect that there are cubic shape-equations for odd p. 

The denominator ^CqSq -l-Pp,-i(co) can be expressed as SQCoifp(co) with Kp 
an even polynomial of order p — 4 in Cq. A recurrence relation for Kp is derived 
in appendix [Cl Apparently Kp{co) is positive for real argument Cq and even 
p > 4. I do not have a proof for general p. However, numerical calculations 
show that up to p = 36 the functions Kp{co) do not have any real zero. 

The solutions of these equations are invariant under rotations. If ai-pm is a 
solution, then also 

a'l;pm = Xl^w("'/3:7)ai;pm' (113) 

m' 

is solution for arbitrary Euler angles a, (3, 7. 



8 Solution of the Shape Equations 

The complete solutions of the shape-equations (jllip is not known to me. After 
all they constitute a set of 2p -f 1 quadratic equations. However, a number of 
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special solutions can be given by restricting to shapes invariant under subgroups 
of the rotation group. Under a subgroup of 0(3) a certain set of functions 
indexed by i 

J2vimCpm{0K)e'""^'' (114) 

m 

is invariant. Thus the amphtude ai-pm is a Hnear combination of the vt^mS, 

i 

The vs should be orthonormal, J2m'^jm''^im = ^ji- Then the shape equations 

read 

^(^iVi,m = (116) 



J2 x/¥TT(-)™- 



77117712 



il,77ll ^12, m2 ' 



P P P 



Multiplication by w*^ summing over m and using orthonormality yields 

aj = '^Wj^i^^i^Ui^ai^, (117) 

jlj2 

Wj,i,,i, = V2P+1 (118) 



77imi77l2 



^ 2^ mo -m M ^ ^^.^^^>..rr>.^^.. 



mi 7722 — m 



If we require Vj ,^ = {—)"^Vj^-m, which implies real a, then we obtain after 
replacing — m by m under the sum 



mmim2 



\/2p+l ( mi ma m ) ^J.™^«i.™i^*2,m2- (US) 



Since the 3j Wigner symbol is invariant under change of the sign of all ms and 
the product of the vs changes into its conjugate complex under this change of 
sign, w is real. Moreover it is invariant under permutation of the indices j, ii, 

12- 

We consider two classes of subgroups of the 0(3). The first class contains 
shapes with an n fold rotation axis, the second one shapes with tetrahedral, 
octahedral and icosahedral symmetry. Obviously the shapes of the last class 
are special shapes of the first class of groups. Without loss of generality I put 
7 = 1 in this section. Otherwise one should read a/7 instead of a. 



8.1 Solutions with rotational symmetry ciround the z-eixis 

If the z-axis is an n fold rotation axis, then the shapes are invariant under the 
cyclic group C„. Since p is even, the shapes have inversion symmetry. Thus 
they are invariant under the groups Cnh for even n and under S2n for odd n. 

Only amplitudes ai-p^kn with integer k running from —[p/n] to [p/n] con- 
tribute. If in particular n > p, then only ai-p^ contributes and since (pO;pO,pO) 
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Table 1 Solutions of the shape equations 
for the dihedral groups Dnh n even 







^l;p,0 




^l;p,2n 


^l;p,3n 




(4 oo) 




2 48576 










(4: 4) 




1 59799 


-0 Q54Q8 








(4 2) 


(2 1 


-0 SQQRO 


-1 26332 


0.71624 




Ox, 


(4 2) 


(2 1 


VJ . -J '.J ^ J. u 


-0 Q8258 


1 29984 






(6 oo) 


(1 0) 


-2 98035 










(6 6) 


(1 1) 


2.70941 


-3 74359 








(6 4) 


fl Tl 


7 45088 


-13 9393 








(6 2) 


(2 2) 


2.70941 


-3.47040 


00000 


-1 40885 




(e, 2) 


(2 2) 


-12.1077 


-4.77180 


8.71207 


-7.07772 




(6 2) 


(2 2) 


-4.40279 


-1 73520 


-1 90082 


1.16987 




(6 2) 


(2 2) 


93136 


-0 95436 


1.04545 


-1.41554 




(6 2) 


(2 2) 


33868 


-0.77601 


-0 63361 


52318 






fl vi\ 


3 40498 










(8 6) 


fl Tl 


1 83345 


-1 24558 








4) 


(2 1 


-0 22918 


-1 48526 


n 99999 






(8 4) 


(2 1) 


2 97936 


-1.12039 


1.70705 






(10. oo) 




-3 78269 










fin inl 


(1 11 

y^T^/ 


7.74691 


-13 5249 








Cm 81 




1 37562 


-1 60639 










n 11 

y'-T'-) 


-1 82513 


-1 31296 








(M) 4) 


(2 1 1 


-139.718 


-123.111 


-48 3395 






flO 4) 


(2 11 


99612 


-1 00375 


-1.19470 






(12.oo) 


(1,0) 


4.12619 










(12,10) 


(la) 


1.34073 


-1.10158 








(12,8) 


(1,1) 


-1.74997 


-2.08835 








(12,6) 


(2,1) 


0.91262 


-1.15888 


0.62268 






(12,4) 


(2,2) 


-5.68198 


-8.02888 


-0.72666 


2.86225 




(12,4) 


(2,2) 


5.08372 


-2.41363 


-1.98403 


5.50958 




(12,4) 


(2,2) 


5.05473 


-1.95363 


1.30127 


-2.91433 


Oh 


(12,4) 


(2,2) 


-0.55198 


-0.93569 


0.89981 


-0.56841 


Oh 


(12,4) 


(2,2) 


-1.64037 


-0.45582 


-1.99424 


0.56834 


Oh 



docs not vanish for even p, there exist for aU p > 2 shapes invariant nnder any 
rotations around the z-axis, that is under Caoh- <p, then several amplitudes 
contribute. Real r implies the condition a\-p-m = (— )™ai;p.m- Distinguishing 
between real and imaginary part of the amplitudes one obtains 2[p/n] + 1 equa- 
tions for 2[p/n] -|- 1 unknowns. By an appropriate rotation around the z-axis 
0'i-p,n can be made real. Then only 2[p/n] unknowns are left. I did not find 
solutions for real ai;p„ 7^ and complex ai-p^kn for 7^ 0,1, but I cannot 
exclude the existence of such solutions. 

If one requires that all a\-p^kn are real, then it is easy to find solutions. They 
are invariant under dihedral groups Dnh for even n, and under Dnd for odd n. 
Shapes were determined up to = 12 and v := \p/n] -|- 1 < 4. 

The number v of independent amplitudes is separated into Vg amplitudes 
o-i;p.kn with cvcu k and i>o amplitudes with odd k. The shape cq. expresses 
the amplitude for even k on the r.h.s. by products of two amplitudes with even 
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Table 2 Solutions of the shape equations 
for the dihedral groups Dnd n odd 



















(1 1) 


-1 06533 


-1.27331 










fl 11 


-1 08376 


-0 86455 










(2 1 1 


n 60209 


-0.91971 


-0 55461 






ffi 31 


(2.1) 


1 3 2460 

A.KJ • ^^xjyj 


-7 99805 


8 38842 




Ox 




(1.1) 


-1 30961 


-1.54101 








f8 51 


(1.1) 


1 83345 


-1.24558 








fS 31 


f2 1 1 


-1 01858 


-1 55591 


54095 






fS 31 


^2 11 


88277 


-1 81826 


2 21627 






fin Qi 


f1 1 1 


-1.72153 


-1.27078 








(10,7) 


(1,1) 


6.00780 


-9.66535 








(10,5) 


(2,1) 


-0.84286 


-1.34290 


-0.73338 






(10,3) 


(2,2) 


2.72446 


-7.51623 


3.45336 


-5.79075 




(10,3) 


(2,2) 


-1.30071 


-0.41447 


-0.92679 


0.90541 


/h 


(10,3) 


(2,2) 


-1.57411 


-0.19827 


0.77110 


-1.02866 


Oh 


(12,11) 


(1,1) 


-2.37941 


-2.98770 








(12,9) 


(1,1) 


-4.89423 


-7.23640 








(12,5) 


(2,1) 


3.54583 


-1.17014 


-0.25936 






(12,5) 


(2,1) 


1.18220 


-0.61091 


0.98335 







fcs and products of two amplitudes with odd fcs, but no mixed contributions, 
whereas the eqs. for the amplitudes with odd ks contain on the r.h.s. mixed 
products of one amplitude with even and one with odd A;. As a consequence to 
any solution with non- vanishing odd amplitudes there exists a second one with 
reversed sign of the odd amplitudes. In the tables only one of these solutions is 
given. 

8.2 Tetrahedral, octahedral and icosahedral groups 

Other subgroups of the 0(3) are the tetrahedral, the octahedral and the icosa- 
hedral groups. Due to inversion symmetry the corresponding groups arc T;j, Oh, 
and Ih- Only floating bodies of equilibrium with mirror symmetry have been 
found, since all these groups have mirror planes. I do not know whether there 
exist any shapes without this symmetry 

Solutions of the shape equations are listed in the following tables. For the 
tetrahedral group one can again distinguish between even and odd contributions. 
The numbers v{p) of linearly independent functions Tk obey 

Voct{p) = t'teteb) = Z^tetob + e), (120) 
Z/oct(p + 12) = !^oct(p) + 1, Z^ico(p + 30) = !/ico(p) + 1- (121) 

The spherical harmonics invariant under these groups are constructed by 

requiring that they are invariant under two elements generating the group. The 
first group element is a rotation about the z-axis by tt for T^, tt/2 for Oh, and 
27r/5 for Ih- The second is a rotation about another axis by 27r/3, 7r/2, and 
27r/5 for Th, Oh, and Ih, resp., described by the rotation matrix D{a,j3,^), 
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Table 3 Solutions of the shape 


equations 






for the tetrahedral group Th 




r) 






n 1 o 




n 1 c 
^l;p,o 










n-i 11-1 
"'l;p,10 


n-i 1 o 
















^l;p,22 










6 


(1 1) 












-^h 


10 




ij . O^O -L -L 


-0 66870 


-0 82Q42 


-0.13114 


\J. t/O 1 ZfVJ 








0.48044 












12 


(2 1] 


-6.12511 


-4 20630 


9 82449 


1 6 5288 


-19 1622 








-10 4288 


n 69076 










12 


(2 1) 




-0 2381 7 


36961 


93590 


14896 










498 '55 








Jx. 


16 


(2 1) 




-0 4881 1 

~\J.*±00 J- J- 


-2.41712 


1 26394 


QQQQ1 








4Q220 


2.10476 


-1 00183 


-1.20704 




Jv. 


16 


(2 1) 


n Q818fi 


-0 40466 


07028 


1.04785 


40520 










54821 


-0 83056 


44283 






18 


(2 2) 


n oo'^fio 


-0 68341 


-0 14607 


-0 6531 8 


61 488 








\J. I \J'-t\J ± 


-0 70238 


-0 46433 


20888 


59331 




18 


(2 2) 


2 S2221 


-3 88457 


-1.49720 


-1 72639 


-0 96421 








-0 nnQQ8 


-2 1 1 050 


-0 1 5926 


-1.76012 


3 00022 

(J . VJ VJ VJ 




18 


(2 2) 




72989 


-0 1 3223 


-0 01 543 

VJ . Vj -L (JtiO 


-0 89094 








68875 


37789 


-0 39434 


-0 58803 


-0 50005 


-ih 


18 


(2 2) 


-4 7'?748 

4:. ( O ( ■rtO 


38021 


2 98822 


-4 42878 


2 25081 








Q 2Q448 


3 98375 


-5.72487 


3 68818 


56796 

VJ . fJVJ 1 tJ\J 




20 


(2 1] 


1 "^8678 


-0.45717 


-0 3551 9 


f) 89628 


f) 49923 










86536 


08246 


1 02995 


-0 82637 








n 40742 










Jx. 


22 


(2 2) 






1 1 5 358 

J_ XfJ . 0<J0 


-273 782 


161.314 








1 04 24fi 


1 28 28^ 


83 0346 


118 399 


-1 35 98^ 








197.709 


Til 876 










22 


(2 2) 


-0.16701 


-1 08554 


-0 64291 

VJ . VJtiZj iJ J_ 


-0 95793 


2 35878 








'?6474 


-0 42041 


29053 


-2 02467 


-0 47579 








1.12764 


9^476 










22 






-9 21 867 


-1 96687 


-4 21 fl'^X 


54457 








-1.31708 


-1 88937 


-0 28461 


-3 3991 1 

O . (J t? t/ -L -L 


30053 

VJ . \J\J\JKJKJ 








-1.11601 


7 22239 










22 


f2 21 


2.63933 


-3.50688 


-1.96687 


-7.01919 


0.54457 








5.59288 


-1.88937 


3.69040 


-3.39911 


-5.87815 








-1.11601 


3.97000 










22 


(2,2) 


0.67761 


0.69725 


-0.19515 


-0.34288 


-0.82641 








0.84351 


-0.26707 


0.48523 


0.03159 


-0.75424 








-0.71327 


-0.39702 










22 


(2,2) 


-1.68119 


2.41122 


1.93518 


-1.18573 


-2.47489 








2.91702 


1.68361 


1.67804 


4.15670 


-2.60831 








-0.22901 


-1.37298 











tetrahedral L»(f,f,0), 

octahedral £)(0, f,0), (122) 
icosahedral Z)(0,/3, |), cos/J=-7=, sin/?=-^. 
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Table 3 cont. 


Solutions of the s 


hape equations 






for the tetrahedral group Th cont. 




D 


V*^ e 7 o J 




fin o 


n ^ A 




rt n o 








n-i 11-1 
"'l;p,10 




rti 1 A 


rt 1 1 (2 


rt 1 no 










^l;p,22 


'^l;p,24 








26 


(2,2) 




-17.3122 














1 Q4fi8fi 

J- • t7*±UOU 


fi Q0372 


5 42601 


4 24205 


1 16469 










-6 975^1 


11.0403 


1 5 8632 






26 




0.44976 


3.76044 


-3.19985 


1.10830 


7.15029 








1.21251 


1.12982 


0.86551 


-4.45935 


0.16510 








-6.08814 


-0.94012 


3.61602 


-2.81306 






26 


(2,2) 


0.88059 


-5.99324 


-2.81215 


-2.56502 


5.52238 








-1.39801 


0.66375 


-0.71139 


-3.85344 


-0.12650 








-5.19225 


0.69595 




4.69009 






26 


(2,2) 


0.89229 


0.89158 


-0.39047 


-0.59796 


-0.44156 








0.86347 


-0.43012 


0.92514 


-0.43095 


0.18639 








-0.46987 


-1.08764 


-0.62933 


-0.44414 







Table 4: Equal shapes 

SPiA = '^£'2,4,1 = '^D3,4 = Sp^j. 

So^fi = SD2,eA^ 

Sdi.c^ = SooXi.! = SD2.6.:i = So.a- 
Sv2A).r> = = Si).jXiA = St.c, = Sjxi- 

SDi^S,2 — Sd^^s,2 — Sq^s- 

1504, 10, 2 = »5d3,10,3 = »5o,10) 
^D5,10 = »^g3,10,2 = St,10 = ^1,10- 
5^4,12,3 = 50,12,2: 
^04,12,5 = So,12,3, 
>5d4,12,4 = <5o,12,l, 
'^D5,12,4 = '^r,12,2 = '^/,12- 



If for Th the amplitudes ai;p,2fc for odd k vanish, then the solutions belong 
to 0;i. I did not tabulate them in the table for the tetrahedral group. If one 
wishes to insert them, then one should multiply the amplitudes a4k by (—)*', 
since different orientations of the equivalent mirror-planes have been chosen for 
Th and Oh ■ For Oh and Ih one can no longer distinguish between even and odd 
amplitudes. 

The solutions which correspond to the groups Oh and Ih arc indicated in the 
tables for the dihedral groups. The solutions invariant under Ih arc indicated in 
the table for the tetrahedral group. Some solutions appear several times in the 
tables in different orientations. A list of equal shapes is givc^n in tabic 4. The 
shapes are denoted by SG,p,i, where G denotes one of the groups I?„, T, O and 
/, resp. and i numbers the solutions for given G and p {i is omitted, if only one 
solution for G and p is listed) . 

Maple was able to solve the set of shape equations without giving any further 
hint apart from the case {ug, t^o) = (2, 2) for the tetrahedral group. Therefore in 
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Table 5 Solutions of the shape equations 
for the octahedral group Oh 



p 


J/ 


ai;p,20 


ai;p,4 
Cll;p,24 


ai;p,8 


ai;p,i2 


ai;p,16 


4 


1 


1.59799 


0.95498 








6 


1 


7.45088 


-13.9393 








8 


1 


2.97936 


1.12039 


1.70705 






10 


1 


0.99612 


-1.00375 


-1.19470 






12 


2 


-0.55198 


-0.93569 


0.89981 


-0.56841 




12 


2 


5.05473 


1.95363 


1.30127 


2.91433 




12 


2 


-1.64037 


0.45582 


-1.99424 


-0.56834 




14 


1 


1.27193 


-0.93532 


-1.00368 


-1.21871 




16 


2 


-1.25843 


-2.46818 


0.71252 


1.71627 


-1.52220 


16 


2 


6.23060 


1.98796 


1.77254 


1.91038 


3.42907 


16 


2 


-7.22633 


3.60347 


-5.11674 


-8.22622 


-1.60500 


18 


2 


-20.7093 


32.0293 


-71.4902 


109.660 


-11.7988 


18 


2 


-3.80429 


2.89479 


-0.31593 


5.60740 


2.23271 


18 


2 


1.04507 


/A A A r\^f\ 

-0.44079 


-f A £^ c\ r\ A 

-1.43304 


0.18341 


-1.13512 


20 


2 


1.59394 


-0.85772 


0.65038 


1.28437 


1.54183 






0.26933 










22 


2 


-1.20555 


0.85156 


-0.10269 


0.83004 


1.41591 






0.57425 










22 


2 


-0.22188 


-0.74186 


2.78358 


-0.47953 


-2.36217 






1.34346 










22 


2 


2.68380 


-0.84569 


-3.04632 


-1.10895 


-0.08719 






-2.72485 










24 


3 


1.18937 


0.54116 


-0.11816 


1.28275 


-0.82343 






0.61481 


0.68382 








24 


3 


1.67166 


0.04356 


r\ -A C\ A ^ A 

0.12474 


1.32901 


0.00398 






1.05199 


0.65127 








24 


3 


1.35890 


-0.55989 


0.33293 


0.71232 


0.94372 






1.02056 


0.27116 








24 


3 


1.85881 


0.00554 


0.86643 


-0.36989 


1.76993 






U.OU ( Uo 


n 770/1 ^ 








24 


3 


0.75146 


0.64693 


0.61258 


-1.06502 


0.92200 






-0.46093 


0.64810 








26 


2 


-1.87265 


1.19575 


0.00293 


0.73397 


1.43595 






1.71925 


0.82798 








26 


2 


-0.08100 


-3.16660 


7.90156 


1.47491 


-4.48413 






-6.19642 


4.25103 








26 


2 


4.63643 


-0.57937 


-5.85327 


-2.88496 


-0.19160 






0.38293 


-5.16865 









this case the equations were rearranged as shown in appendix[Bj Apparently the 
case (2, 2) is easier solved for the dihedral groups, since some of the coefficients 
"^jMM vanish in this later case. 

In general 2"^ solutions were obtained including the trivial one, for which 
all amplitudes vanish. Only real solutions are listed, since in the calculation 
ai;p,-m was expressed by (— )'"ai;p,m, whereas it should be (— )™a*.p 
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Table 5 cont. Solutions of the shape equations 
for the octahedral group Oh 





1/ 


'^l;p,0 




^l;p,8 


^l;p,12 


16 










^l;p,28 






28 


3 


-11.2790 


-15 3654 


-2 10025 


-3 08484 


3 01089 

KJ -yj A. yjij ^y 








1 63605 


-11.7977 






28 


3 


-8 1Q3fi2 


-77.0755 


-6 62562 


9 30006 


32 7654 






48 2929 

TTl^ * \J \J 


35 2893 


-41 3066 






28 


3 


1.31496 


-0 12505 


-0 16182 

Vy ■ j_ vy j_ Ly ^ 


1.13151 


54722 






-0 02126 


1.04212 


39059 






28 


3 


4.96143 


2 931 82 


3 87509 


-3 45602 


33561 






4 06348 


-1.77305 


3 70403 






28 


3 




2 16990 




-1 80224 


1.82715 






5 30833 


32544 


4.48484 






30 


3 


-10 3305 


3 75338 


5.18714 


4 88790 


4 09221 






4,37127 


4 621 63 


7 22909 






30 


3 


-4 12966 


1 90227 

J- . tJ\j i-ii-i 1 


1.61647 


-0 22958 


5 291 45 






-0 70240 


3 63614 


2 20732 






30 


3 


-0 32893 


-0.47577 


1 23931 

X • ^ O ly O -i- 


1.34217 


-1.44945 






41 590 


-1.13186 


1.04467 






30 


3 


80723 


36595 


-1.83491 


-0.45764 


-0 88932 

W . t/ ry ij 






1 37893 


22661 


-1 34800 






30 


3 


73.5784 


-105.971 


206 495 


-395 192 


731 230 






-842 698 


143 693 


7.16578 






32 


3 


-0 94505 


-1 62873 


1 80496 


-2 69679 


-1 24200 






1.76241 


2 78195 


-1.70814 


-0 86640 

V.y . (_J Vy t: W 




32 


3 


2 02108 


-0 08638 


-0 29673 


1 28426 


1.07854 






37970 


1 9005 


1.46548 


62050 




32 


3 


-1 16350 


1 20893 


-0.44735 


08591 


-0 46335 






-1 25548 


-1.57057 


-0 62832 


0.12415 




32 


3 


7 53901 


1 2 2486 


8 93730 


-5 85501 


-5 29201 

^ m ^ ^ A. 






-0 18349 


1.41254 


-7.58121 


1 0563 




32 


3 


47 6029 


30 0645 


24 1029 


-3 30451 


-1.96647 






7 99337 


1 1 6529 


-3 80089 


35 2026 




34 


3 


-11.4914 


-9 42515 


25 0654 


30 3572 


-9.41248 






-9.23599 


1.92303 


-17.1336 


25.2407 




34 


3 


4.10996 


-1.74467 


-1.65191 


0.50294 


-3.74252 






-2.26974 


0.11285 


-3.55783 


-2.00760 




34 


3 


-5.79444 


0.80954 


4.16799 


5.11033 


0.56450 






0.94328 


2.59973 


0.35243 


5.36668 




34 


3 


0.51337 


0.24287 


-1.08014 


-0.06906 


-0.82773 






0.02861 


0.97644 


-0.20870 


-0.77065 




34 


3 


1.09993 


-0.90039 


0.44903 


-0.02780 


-0.34114 






-0.84004 


-1.20752 


-0.96928 


-0.08435 





Case 1^ — 2 As an example the set of equations (1117p for the case = 2 is 
considered, 

ai = wiiial + 2wii2aia2 + wi22al, (123) 





Table 6 Solutions of the shap 


e equations 






for the icosahedral group Ih 




r) 


y 


t*l;p,0 


n -1 r 


n 1 in 

'^l;p,lU 




(7,1 on 






""lip, 25 


f^l;p,30 


"'l;p,35 






6 


1 


-1 08376 










10 


1 


-f) 84286 


-1 342Q0 








12 


1 


1 18220 


-0 filOQI 


Q8335 






16 


1 




1.87102 


-2 40430 


-1.86451 




18 


1 


-1.24551 


0.70244 


-0 56268 


1.11003 




20 


1 


3581 fi 


1.06413 


1.72730 


53368 


0.94937 


22 


1 


-9.42141 


-3 84Q03 


6 02057 


-6 60823 


-6.45637 


24 


1 


1 28346 


-0 8441 5 


n 57494 


-0 5Q346 


1.25220 


26 


1 


-3 3536Q 


-6 46073 


-4 1 7676 


7 321 38 


3 70065 






5 20342 










28 


1 


-21 27Q0 


-4 33343 


10.4040 


-11.4419 


12.2013 






14.1643 










30 


2 


-1 2Q745 


1 00742 


-0 67877 


n 53443 


-0.66788 








-0 00751 








30 


2 


435Q7 


-0 78Q87 


-0 Q0073 


-1 69849 


-0.33963 






-0.64517 


-0 90325 








30 


2 


3 Q4227 


-1.15076 


6 83Q83 


4 80458 


4.41657 






-3 58280 


3 85631 








32 


1 


-25 3074 


-34 3654 


-6 83507 


30 6595 


-30 5265 






-22.3716 


-28.6144 








34 


1 


7.10336 


0.32668 


-2.79493 


3.20114 


-3.38581 






3.53378 


4.72435 








36 


2 


-l.U ( ^Oi 


z.uoo / u 


1 . ( oizy 


1 .UOOO 1 


-3.30158 






-0.87907 


-2.15745 


-2.00621 






ou 


2 


-0.83679 


-0.11942 


-1.82417 


-0.21073 


1.13536 






1.30552 


-0.61666 


1.06790 






36 


2 


1.28515 


-1.18256 


0.83404 


-0.59249 


0.54385 






-0.77680 


1.60342 


-0.01678 










a-i = 




2wi22Q:ia2 - 







The ratio 



a2 



is introduced for the non-trivial solutions. Then the equation 

W\i2 + 2wi22r + W222r'^ 

r = 

Will + 2tUii2r -I- Wi22?'^ 

is obtained, which rewrites as an equation of third order in r, 

wi22r^ + (2wii2 - W222)r^ + {win - 2^122)^ - W112 = 0. 
Thus there are three non-trivial solutions. One or three of them are real. 



(124) 
(125) 

(126) 
(127) 
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9 Higher Orders in the deformation 



As in the preceeding section only central symmetric bodies are considered. Then 
the expansion of r and h in spherical harmonics contains only those with even 
I. In higher orders the equation for Van can be used to determine hn{P,a) up 
to an /3, a- independent contribution hn-oo- The difference — CoroVan yields 
the condition 

(128) 

for I > 0, where In-im contains only terms a™ with m < n — 1. For ^ 7^ p it 
yields 

an,im = — — {In;lm — 2^3; (aia„_i ) ;(,„ ) . (129) 
All 

For / = p one obtains a set of linear eqs. for the a„_i.pm 

^n—l;pm 2 ^ ^ ^^mm' 0'n — l\pm' — ^n;m (1^^) 



Wl 



th 



Mjnm' = {pm]pm- m' ,p'm') °'^'^"^ " , (131) 

7 

J ^n — 1:00 -^n;pm 



^1;00 So'll;OO^p,-l(co) 

cqSo + 2Pp _i(co) 



(132) 



-r / (j)Q]pQ,lQ){prn]prni,lm2)ai.prnian-i-irn2- 

This set of equations has a unique solution, if the matrix M does not have an 
eigen-value 1/2. An eigen-value 1/2 appears, if the shape-eqs. have a double- 
solution. This means that besides ai-p, also ai.p, + 5ai.p, is solution to first order 
in 5ai;p., which yields the equations 

<5ai;pm — 2 Mmm'^0'l;pm' — 0. (133) 

This is the case, if one does not fix the orientation of the body, since then the 
change of ai;p. by infinitesimal rotations yields such solutions ai;p. -\-5ai-p,. Here 
we have fixed the symmetry axes in the calculations. Thus this case does not 
appear. We cannot exclude that in special cases a double-solution appears. We 
have determined the eigen-values of M for all shapes listed in the tables. No 
eigenvalue 1/2 appeared. The closest one were 0.4541 for 5^4, 12, 5 and 0.5979 
for 5t,26,2- Always one of the eigenvalues equals 1 for the eigen- vector ai;p.. 

9.1 Example for double-solution 

We give an example for a double-solution by further considering the case v — 2. 
Starting from eq. (|127p and denoting the solutions by ri, r2, one obtains 

Will = ^122(2 -I- rir2 + rirs + r2r3), (134) 
W112 = Wi22rir2r3, (135) 
W222 = Wi22(?'i + ?'2 + ^"3 + 2rir2r3). (136) 
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Considering the solution r = ri one obtains 

1 n 

"1 = 77, OL-i 



W122N' W122N' 

N = 2 + rir2+rir3 + r2r3 + 2r^r2r3 + r2. (137) 
This yields the matrix M with the matrix-elements 

Mij = Wyiai + Wjj2a2, (138) 

,f 2 + rir2 + rira + r2r3 + r^r2r3 

= > (139) 

M,2 = M2i = !:i(i±!^, (140) 
M22 = l + ri{ri+r2+T:i) + 2rlr2r:, _ ^^^^^ 

The eigenvalues of this matrix are 

1 + rir2 + rir3 + r?r2r3 
A*i = 1, M2 = • (142) 

One finds 

Thus one obtains an eigenvalue /i = 1/2 if and only if ri equals one of the other 
solutions r, . 



9.2 Reparametrization 

By iterating the expansion in e it became apparent that the coefhcients r„:oo 
and hnfio were undetermined, whereas all other coefficients r„;;m and hn-im were 
determined. 

The reason is the following: One is free to reparametrize r by introducing a 
new r = rsr(e) with a function Sj.(e) that has a Taylor expansion in e. Simul- 
taneously ft. has to be changed to h — hsr{e). Thus one is free to choose r„.oo 
arbitrarily. One may for example choose r„:oo — for n > or require a volume 
independent of e. 

If the r„.ooS are fixed one is still free to reparametrize h to h = hsh{e), where 
again Sft(e) must be Taylor expandable. A possible choice is hn-QO = for n > 1. 



Acknowledgment I am indebted to Christian Wegner for useful discussions 
on section 3. 



A Supplement to section [3] 

A.l Ultraspherical (Gegenbauer) polynomials 

For more information on ultraspherical polynomials see [181 119j . 
Introducing coordinates 



d 

cos 9, 

k=i+l 



n sin^fc, (144) 
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where 0k runs from to tt for fc > 1 and 6*1 = 0, tt alfows to write the ultrasphe- 
rical harmonics as 

^nd,mci-i,...,m2,mi (^f^' ■■■^l) 

= 'Od / , T^^m.^i ...,m2,mi(^rf-l.--^l): (145) 

/ Ar(md-l+a/2-l) 

rW(^i) = ^cos™^^i (146) 

with the restriction irik > mk-i and mi = 0, 1. For given d and uid these are 
{2md + d— 2){md + d — 3)1/ {{d — 2)!md!) ultraspherical harmonics. They form 
a complete orthonormal set of functions and obey eq. p9|) . where n = md and 
TO stands for the set {nid-i, ...,mi}. 

Gegenbauer's addition theorem [5D] reads 

(n + a)T{a)G'C^ [xx' + ^(1 - •^^)(1 " 2;'2) cos^) 

_ " (l-^2y/2(i_^,2y/2 

- TTSTI) '^n-i (2;)G„_; (.xj 

X + l/2)r(a- 1/2)g|""'/'^(cos(/.), (147) 
with the norm 

da;(l - x2)"-i/2G(f ) (a;)G(;^) (x) = 7V("),5„,„, (148) 

n!(n + aji ^(a) 



Repeated use of this addition theorem with u = xed + \/\ — x^u', unit vector, 
X = cos 6 d{vL), similarly for v, and u' • v' = cos0 yields eq. (fT8|l . 



A. 2 Determination of the expansion coefficients c„ 

The Legendre polynomials are ultraspherical polynomials, Pn{x) ~ Gn^'^\x). 
Thus it is sufhcient to consider the expansion eq. (fT7|) . Orthogonality of the 
polynomials p48p yields 

CrMT^'^ = 2 f' dxx{l-xy^'-'^'GiT''\x). (150) 
G is represented by Rodriguez' formula 



Ad/2~i) _ ^ \ — 



T{^)T{2n + d-2) 
22" (2n)!r(d- 2)r(2n+ ^) ' 



(152) 

Thus 



"1 ^2n 
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Partial integration yields [T3] 



i(<i/2+l 

'^2n-2 

With 



ffe)Grrr^(o). (154) 

2n-2 



G^n^^(0)^("+^'/'-') (155) 



one obtains the expression for c„ in eg. p7|) . 

B The case (z/c, Vo) = (2, 2) 

The shape-equations can be written 

= c^~c'^^^, (156) 

= Cy~c'y^^, (157) 

= CuuU + CuvV, (158) 

= CuvU + CmV. (159) 

with 

= a; - woooa;^ - 2woo2Xy - wo222/^, (160) 

Cy ^ y ~ woQ2x'^ - 2wo22xy - ^2222/^, (161) 

4™ = WQiiu'^ + 2woi3UV + W033v'^ , (162) 

= W2iiu^ + 2^213™ + ^«233^'^, (163) 

Cuu = 1 - 2woiia; - 2w2ii2;, (164) 

Cuv = -2woi3x - 2w2i3y, (165) 

Ct,„ = 1 - 2wo33a; - 2w233.y- (166) 



Eqs. USHIIIMI) yield 



^2 



cL = 0., (167) 



V 

Insertion in (|162|163p yields 



(168) 



uv 

Cxuv, Cxuv = WQiiCyy " 2u'oi3C„i, + ^033 C«« , (169) 



I uv 

Cyuv = Cyuv, Cyuv = W2llCvv - 2w21?,Cuv + W233C««. (170) 

and substitution in eqs. ()156|157p and elimination of ^ yields 

0. (171) 

The eqs. p671 117ip allow the determination of x and y. 
From (|156|169|168p one obtains 

C-x^vv ^xuv^ — (1'^^) 

which allows the calculation of u. Finally eq. (|158p yields v. One notes that 
the solutions with u ^ v ^ have been eliminated by use of (|1671 1168P . They 
can be obtainded from Cx = Cy — 0. 
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C Denominator of the coefficient 7 

The Legendre functions obeys the recursion relation 

(p + 2)Pp+i,_i(a:) - {2p+l)xPp,^i{x) + (p- l)Pp_i,_i(x) - 0. (173) 

Using it also forp ~f p±l, one can eliminate Pp±i^-i{x) and obtain the recursion 
relation 

(p + 2)(p + 3)(2p-l)Pp+2,-i(x) 
+ (2p + 1) (2p2 + 2p - 3 - (2p - 1) (2j9 + 3)x^)Pp^-iix) 
+ (p-l)(p-2)(2p + 3)Pp_2.-i(a;)=0. (174) 



Substituting 



isoco + Pp _i(co) = soCo(l - cg)i^p(co) (175) 



one obtains with Co = a; 



2 

One easily calculates 



{p + 2)(p + 3){2p-l)Kp+2{x) 
+ {2p + l)(2p2 + 2p - 3 - (2p - l){2p + 3)x^)Kp{x) 
+ ip~l)ip-2)i2p + 3)Kp^2{x) 

\2p- l)(2p+l)(2p + 3). (176) 



K2=0, K4 = l- (177) 



The recursion relation (|176p yields 

Ke = ^(l + Ux^), 
16 

Ks = ^(9-26a;2 + 65a;4), 

Kio = 7^ (5 + 89x^ - 289a;'' + 323a;^) , 
256 

K12 = Y^(''43- 5032a;2 + 34238a;^- 72353a;^ + 52003a:*^). (178) 

Numerical calculations yield no real zeroes of Kp{x) up to p = 36. Probably all 
Kp(x) are strictly positive for real x. 

References 

[1] R.D. Mauldin (ed.), The Scottish Book, Birkhiiuser Boston 1981 

[2] H. Auerbach, Sur un probleme de M. Ulam concernant I'equilibre des corps 
flottant, Studia Math. 7 (1938) 121-142 

[3] F. Wegner, Floating Bodies of Equilibrium Studies of Applied mathematics 
111 (2003) 167-183 



28 



[4] F. Wegner, Floating Bodies of Equilibrium. Explicit Solution, e-Print archive 
|physics/0603160| 

[5] F. Wegner, Floating Bodies of Equilibrium in 2D, the Tire Track Prob- 
lem and Electrons in an Inhomogeneous Magnetic Field, e-print archive 
^physics/ 0701241 

[6] P. Bouguer, Traite du navire, de sa construction et de ses mouvements 
Jombert Paris (1746) 

[7] C. Dupin, Applications de Geometric et de Mechanique, a la marine, aux 
ponts et chaussees, etc Bacheher Paris (1822) 

[8] P. Appell, Traite de Mecanique Rationelle, Gauthiers-Villars, Paris (1921), 
tome 3 (1921) 194-233 

[9] A.G. Webster, The Dynamics of Particles and of Rigid, Elastic, and Fluid 
Bodies 2nd ed., G.E. Stechert, New York (1922) 

[10] E.N. Gilbert, How Things Float, Amer. Math. Monthly 98 (1991) 201-216 

[11] R. Schneider, Uber eine Integralgleichung in der Theorie konvexer Korper, 
Math. Nachr. 44 (1970) 55-75 

[12] R. Schneider, Functional Equations Connected with Rotations and Their 
Geometrical Applications, L'Enseignement Math. 16 (1970) 297-305 

[13] K.J. Falconer, Applications of a Result on Spherical Integration to the The- 
ory of Convex Sets, Amer. Math. Monthly 90 (1983) 690-693 

[14] E.P. Wigner, Gruppentheorie and ihre Anwendung auf die Quanten- 
mechanik der Atomspektren, Vieweg, Braunschweig, Germany, (1931) 

[15] E.P. Wigner, Group Theory and Rs Application to the Quantum Mechanics 
of Atomic Spectra, Academic Press, New York, (1959) 

[16] A.R. Edmonds, Angular Momentum in Quantum Mechanics, Princeton 
University Press (1957) 

[17] L. Montejano, On a problem of Ulam concerning a characterization of the 
sphere. Studies in Appl. Math. LIII (1974) 243-248 

[18] M. Abramowitz, LA. Stegun (eds.). Handbook of Mathematical Functions, 
chapter 22, Dover 

[19] A. Erdelyi (ed.), Bateman Manuscript Project: Higher Transcendental 
Functions, volume I, section 3.15, McGraw-Hill 

[20] L. Gegenbauer, Das Additionstheorem der Funktionen C^(a;), Sitzungsber. 
der math.-nat. Kl. Akad. Wiss. Wien Ha, 102 (1893) 942-950 



29 



